Nonclassicality of a photon-subtracted Gaussian field 
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We investigate the nonclassicality of a photon-subtracted Gaussian field, which was produced in 
a recent experiment, using negativity of the Wigner function and the non-existence of well-behaved 
positive P function. We obtain the condition to see negativity of the Wigner function for the 
case including the mixed Gaussian incoming field, the threshold photodetection and the inefficient 
homodyne measurement. We show how similar the photon-subtracted state is to a superposition of 
coherent states. 
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I. INTRODUCTION 

The recent development of quantum optics has opened 
the possibility to generate and manipulate various non- 
classical light fields, which cannot be described by clas- 
sical theory, in a real laboratory. It is generally accepted 
that the presence of a positive well-defined P function (a 
quasiprobability function in phase space |3|) signals the 
field classical Q ; otherwise the field is categorized as non- 
classical. A stronger constraint on nonclassicality is the 
presence of negativity in the Wigner function (another 
quasiprobability function) of the field [f| ■ While a Gaus- 
sian field may not have its P function, its Wigner func- 
tion never becomes negative. For example, the squeezed 
vacuum state is represented by its Gaussian Wigner func- 
tion while its P function does not exist [fj. It is also 
known that a Gaussian field remains Gaussian by lin- 
ear transformations which correspond to basic tools in 
a quantum optics laboratory such as a phase shifter, a 
beam splitter and a squeezer 0, ■ 

Two better-known nonclassical fields are a squeezed 
state and a superposition of two separate coherent states 
(coherent-state superposition). The two kinds of states 
are closely related to probably the most fundamen- 
tal and intriguing paradoxes in quantum theory, i.e., 
the Einstein-Podolsky-Rosen paradox Jjj for a two-mode 
squeezed state and the Schrodinger's cat paradox 8] for 
a coherent-state superposition. They are also known as 
useful resources for various schemes in quantum informa- 
tion processing. A squeezed state and a coherent-state 
superposition manifest different types of nonclassicality. 
Whereas a squeezed state is a Gaussian field, a coherent- 
state superposition is non-Gaussian and shows a large 
amount of negativity in its Wigner function. There was 
an early attempt to relate the two states through quan- 
tum noise of arbitrary strength 0. Dakna et al. [To| 
considered a connection between the two states by sub- 
tracting a precise number of photons from a squeezed 
field. They also showed that any quantum state can be 
generated from the vacuum by application of the coher- 
ent displacement operator and adding photons 0. On 
the other hand, it has been reported that by squeezing 



a single photon state one can generate a state which has 
almost unit fidelity to a coherent-state superposition of 
small amplitude [l2 ]. 

It is only very recently that a traveling non-Gaussian 
field was experimentally generated by subtracting a pho- 
ton from a squeezed vacuum by Wenger et al. ll.'Sl ) . They 
used a beam splitter and a threshold detector to sub- 
tract a photon from the squeezed field, but the recon- 
structed Wigner function failed to show a negative value 
[lflj . It is, thus, timely to analyze the generation of a 
non-Gaussian state in relation to the status of experi- 
ments. In particular, as such the state forms a start- 
ing point for distillation of a continuous- variable field for 
quantum information processing |14| and may improve 
the efficiency of quantum teleportation the study 
will be of use. In this paper, we assess the nonclassical- 
ity of a photon-subtracted Gaussian field and study how 
similar this state is to a coherent-state superposition. We 
assess negativity of the Wigner function in conjunction 
with the non-existence of the positive P function. 

II. FIELD GENERATED BY SUBTRACTING A 
PHOTON 

We would like to consider what kind of state one pro- 
duces by eliminating one photon from a simple Gaussian 
function. A single-mode Gaussian field of its density op- 
erator p may be represented by the Weyl characteristic 
function [l6f defined as C(£) = Tr[£)(£)p]: 

C(0=exp(-^ 2 -f^), (1) 

where A and B are determined by the quadrature vari- 
ances of the field. The displacement operator has been 
defined as D(£) — exp(^a t — £*a), where a and a* are 
bosonic annihilation and creation operators respectively. 
Note also that the density operator can be obtained from 
the characteristic function as 

P = ~ / d 2 ZC(0D(-0, (2) 
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which can be straightforwardly obtained using identities 
(1/tt) / d 2 a\a){a\ = 1 and 20] 
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where \a) is a coherent state of amplitude a. Even 
though Eq. (JTJ does not represent a very general Gaus- 
sian field, rotation and/or displacement operation brings 
any Gaussian field to this form. It is useful to start with 
(|T|l because it is extremely challenging to produce a pure 
squeezed state with AB — 1 and the characteristic func- 
tion allows us to treat a single-mode Gaussian state 
of a mixed state. The uncertainty relation is given by 
AB > 1 and the Gaussian state is called squeezed when 
either A < 1 or B < 1. 

Let us consider the experiment by Wenger et al. |13| . 
First of all, they produce a squeezed Gaussian state then 
this passes through a beam splitter with its transmittiv- 
ity T = t 2 , where the other input port is assumed to be 
served by a vacuum. At the one output of mode 2, we 
conditionally measure a one photon state 1 1)2- The state 
obtained at the other output port of mode 1 was what 
Wenger et al. produced as a non-Gaussian field in their 
experiment. We will evaluate the Wigner function for 
this field of mode 1. 

By beam splitting the squeezed Gaussian field whose 
characteristic function is written as and the vacuum 
of its characteristic function, C„(0 = cxp(— i|£| 2 ), the 
characteristic function for the output field of modes 1 
and 2 is EH 



C<mt 07)0 = exp ( -ixVx 5 



where x = (r) r , rji, £ r , £j) and the correlation matrix 
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V = 
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with 

m=TA+R, n 2 = TB + R, a=tr(A-l), 
c 2 =tr{B-l), m 1 = RA + T, m 2 = RB + T (5) 

and T — t 2 and R = r 2 . 

We then use the two-mode version of f° r the density 
operator of the output field: 



Pout 



C ou t(vA)D 1 (-v)D2(~0d 2 vd 2 ^ (6) 



The density operator for the field of mode 1 conditioned 
on one-photon measurement for mode 2 is 



Pi = N 2<1|/W|1)2. 



(7) 



Throughout the paper, J\f denotes the appropriate nor- 
malization factor. For the case of Eq. |J7J|, 
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[(mi + l)(m 2 + 1)] 3/2 



(8) 



2 (l|Tri[p olit ]|l) 2 2(mim 2 - 1) 

With the knowledge of the one-photon Fock state expec 
tation value of the displacement operator [23, [2^J| 



(l|£(-£)|l)=e-^(l-|£| 2 ), 
the density operator is found to be 



Pi = — 



c(77,0A(-»?)e-— (i-ier)dV 2 e 



The characteristic function is then easily obtained using 
the identity Tr[D(()D(-r])] = ttS^K - rj): 



Ci(C) = 



1 - 



c 2 (m 2 + 1)C 2 
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By Fourier transformation of the Weyl characteristic 
function |2^|, we obtain the Wigner function. Now, the 
first point we are interested in is the negativity of the 
Wigner function. It is clear that the Fourier transform 
of @ has the largest negativity (if any exists) at the ori- 
gin of phase space and the value of the Wigner function 
at the point is 

W ^ X (T + 1)B + R + (T + 1)1 + R (10) 
which has been obtained by substituting the parameters 



(J5J. It is obvious that if A > 1 or B > 1, i.e., the incom- 
ing Gaussian field is not squeezed, W(0) is positive every- 
where. In order to find the exact condition for negativity 
in the Wigner function, we assume that A < 1, B > 1 
and introduce positive parameters x = (A + — A) 

and y = (B + 1)/(B-1). Then the right-hand-side (RHS) 
of HI L)|l becomes 

2T - x + y 
(T-x)(T + y) 
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whose denominator is always negative. The numerator 
becomes positive when the transmittivity satisfies 

AB-l 

r> (l-A)( B -l) <"> 

which always holds when the incoming Gaussian field is 
pure AB = 1 (in other words, if the incoming Gaussian 
field is a pure squeezed state, the Wigner function always 
shows negativity by subtracting a photon from it). 

The P function of the field may be obtained using the 
relation between its characteristic function C{ and the 
Weyl characteristic function [23| : 



CW(C) = C(C)e4 



ICI 2 



(12) 



With use of the characteristic function © and general 
Gaussian integration, we find that the P characteristic 
function is integrable when (n.j — l)(mj + 1) — cf > for 
i = 1,2. By substituting the parameters 10, we find the 
condition equivalent to 2T(A-1) > and 2T{B-1) > 0. 
So if the incoming field is squeezed, it is not possible to 
integrate and no P function exists. Considering the 
positivity of the P function, after a little algebra with 
Fourier transformation of the P characteristic function, 
we find that the P function is positive as far as it exists in 
this case. We conclude that the single-photon subtracted 
field is nonclassical (in the sense of a lack of an accept- 
able P function) provided the original incoming field is 
squeezed. (However, the Wigner function does not nec- 
essarily show negativity for all those nonclassical states 
unless the incoming Gaussian field was pure.) Unless the 
incoming Gaussian field is nonclassical we cannot gen- 
erate a nonclassical state by subtracting a photon from 
it. 

This seemingly trivial result is not obvious at all as 
contrasted by the nonclassicali ty o f a field by adding a 
photon into a Gaussian field (25ll26| . In distinction to the 
case of subtracting a photon, the photon-added Gaussian 
state alw ays shows negativity at the origin of the phase 
space |2(||23,|28|. By adding a photon, a highly classical 
state such as a high-temperature thermal state becomes 
non-classical, showing negativity in its Wigner function. 
The realization of such a photon added state is beyond 
the scope of the current work but we may think of a 
possibility within cavity quantum electrodynamics or the 
phonon state of a driven ion in a cavity |28| . 

We now introduce the coherent-superposition state [24[ 



where the subscripts refer to the respective states. 

The fidelity between p\ and the coherent-state super- 
position 1)13(1 has been calculated from Eqs. JU, $Y.\\. 
and l|14|) and plotted in Fig. The incoming Gaus- 
sian field has been assumed a pure squeezed field. In 
Fig. the solid line is the optimized fidelity between 
the photon-subtracted state and the ideal coherent-state 
superposition by an ideal single-photon detector. The 
fidelity is very high as T > 0.99 regardless of the trans- 
mittivity of the beam splitter when an ideal single-photon 
detector is used. The optimized amplitude of the ideal 
coherent-state superposition is a = 1.16 for the transmit- 
tivity close to unity. If the transmittivity T gets smaller, 
the amplitude of the ideal coherent-state superposition, 
which maximizes the fidelity, also becomes smaller. For 
example, the amplitude will be a = 1.02(1.09) for T = 
0.8 (0.9). However, the fidelity is not sensitive to the 
transmittivity of the beam splitter as shown in Fig. ^ 
because the single photon detector successfully subtracts 
only one photon from the Gaussian state regardless of the 
transmittivity of the beam splitter. In fact, the fidelity 
gets slightly better as the transmittivity becomes smaller, 
due to the fact that both of the states are reduced to the 
exact single photon state as T — > 0. 

It is interesting that the fidelity between the photon 
subtracted field pi and the coherent-state superposition 
is very high. This could have been guessed from their 
photon number distributions. The squeezed vacuum is 
a state with only an even number of photons |6j while 
the coherent-state superposition (|13|l is a state with only 
an odd number of photons [12J . By subtracting one pho- 
ton from the squeezed state, the two states may become 
closer to each other. We see that the photon-subtracted 
squeezed field is close to the coherent-state superposition 
of small amplitudes. One reason can be found again in 
their photon number distributions. The photon number 
distribution of \ip) peaks around \a\ 2 while that of p\ is a 
monotonous decreasing function with regard to the pho- 
ton number. Thus, when a is small, the distributions 
become similar to each other. Of course, this check of 
the photon number distributions gives only a hint as the 
photon number distribution does not necessarily convey 
all the coherence properties of a quantum field. 



m=Af(\a)-\-a)), 



(13) 



where Af = l/yl — exp[— 2a 2 ], to assess its fidelity to 
the photon-subtracted Gaussian state. It is straight- 
forward to calculate the characteristic function of the 
coherent-state superposition from Eq. (jT3|) [|| • The close- 
ness of two states, one of which is a pure state \<p) and 
the other (pure or mixed) is represented by its density 
operator p, is measured by fidelity T: 

T=W\$l =~ fd 2 (C^()C p (0 (14) 



III. EXPERIMENTAL REALITY 

As can be seen in Fig. 1 , the fidelity between the ideal 
coherent-state superposition and the photon-subtracted 
state is not so sensitive to reflectivity of the beam split- 
ter. This seemingly good result is due to an ideal sin- 
gle photon detector assumed for the photon-subtracted 
state pi. As mentioned, the state (J7J) is what is wanted 
to achieve but the available high-efficiency photodetec- 
tor is not able to discern 1 and any number of photons. 
Thus the state experimentally generated using such the 
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FIG. 1: The fidelity between the photon-subtracted state and 
the ideal coherent-state superposition with an ideal single- 
photon detector (solid line) and a threshold detector (dotted 
line). The initial squeezing parameter is exp(2s) = 2.36 and 
the x axis is the transmittivity of the beam splitter T = t . 
The amplitude a of the ideal coherent-state superposition is 
optimized for the maximum fidelity. The optimized amplitude 
a ranges between 1.02 (when T = 0.8) and 1.16 (when T — > 
1). 



threshold photodetector is 

oo 

Pa = N/, 2 (n | pout \n) 2- 



(15) 



Consider the density operator for mode 1 of the output 
field 



Pt = Tr 2 [p out ] = ^ 2{n\p out \n) 2 - 

ra=0 

It is then clear from Eqs.lHSJ and (JTHJ) that 
Pa =M{p t - 2 (0|p out |0) 2 ) 

where 

Pt = \ J C out (v,0)Di(-V)d 2 V 

and 



(16) 



(17) 



(18) 



2 (0|/w|0) 2 = i J C out (r,,Oe-^ 2/2 Di(~v)d 2 vd 2 ^ 

(19) 

Using C(rj, £) we have already discussed, we find the char- 
acteristic function C a (C) f° r Pa- 



C Q (C) = Me-UmC 2 r+n 2 & 



i/(mi + l)(m 2 + l) 



x exp 



bmi+l)^ 



2(mi + l) >r 2(m 2 + l) 
The normalization factor is calculated as 



(20) 



N = 



V(m 1 + l)(m 2 + l) 
A /(mi + l)(m 2 + l)-2' 



The Wigner function obtained by Fourier transformation 
of the characteristic function l|2U|l is what Wenger et al. 
would have reconstructed if the detection efficiency of 
their experiment had been perfect and the modal purity 
unity. 

Let us next consider the negativity of the Wigner func- 
tion. By inspection of the characteristic function, we re- 
alize that the Wigner function has the largest negativity 
(if any) at the origin of the phase space and the value of 
the Wigner function at this point is 

- 2Af f 1 



W a (0) 



7T I ^n x n 2 



| (21) 



i/[ni(mi + 1) - c1}[n 2 (m 2 + 1) - c\ 

By partly substituting the parameters |5"jl. we find that 
the Wigner function becomes negative when 



> 



1 



(22) 



y/{m + A)(n 2 + B) ' V n i n 2 
which becomes a criterion for the transmittivity 

T > 4 -(- 4 + 1 )( j3 + 1 ) (23) 
1 > 3(A-1)(S-1) ' { 6) 

For a pure squeezed Gaussian incoming field, the con- 
dition becomes T > 1/3. It is interesting to note that 
regardless of the degree of squeezing (provided it is not 
zero), we can see the negativity in the Wigner function 
provided the transmittivity is larger than 1/3. 

Let us assess the degree of nonclassicality by the P 
function criterion. With use of the relation (|12f> between 
the characteristic functions, we note that the P charac- 
teristic function for p a is integrable when T(A — 1) > 
and T(B - 1) > 0. We have checked that the P func- 
tion is semi-positive when it exists and conclude that, for 
nonzero transmittivity of the beam splitter, iff the incom- 
ing Gaussian field is squeezed, the any-number photon 
subtracted state p a is nonclassical. Again, the P func- 
tion criterion is weaker than the negativity criterion for 
the Wigner function. 

We now consider how close the field obtained using 
the threshold detector to the coherent-state superposi- 
tion (|13l) . The optimized fidelity has been calculated 
using Eqs. ijTSl Hi)) and QJJ, and plotted in FigUJ It 
tells us that the state which is obtained by subtracting 
any number of photons is similar to the coherent-state 
superposition only when the transmittivity of the beam 
splitter is very high. For example, the fidelity is higher 
than 90% when T > 0.87. In this case, the chance of one 
photon subtraction is more likely. Note that the opti- 
mized amplitude a ranges between 1.02 (when T = 0.8) 
and 1.16 (when T — > 1) in Fig. [J 

A. Inefficient detection and modal purity 

Homodyne detection may be used to reconstruct the 
Wigner function for the field p a . Even though homo- 
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FIG. 2: (Color online) (a) The Wigner function W(z) of 
a photon-subtracted Gaussian state with a threshold detec- 
tor for photon subtraction and ideal homodyne detectors 
for reconstruction of the Wigner function when T — 0.88 
and exp(2s) = 2.36. The minimum negativity is found as 
W(0,0) = -0.52. (b) The Wigner function of a photon sub- 
tracted Gaussian state under the same condition as (a) but 
with homodyne efficiency r\ = 0.75. The minimum negativ- 
ity has been reduced to —0.15. (c) The Wigner function of 
a photon subtracted Gaussian state under the same condi- 
tion as (a) but with homodyne efficiency r\ = 0.75 and with 
the modal purity factor 0.7. The negativity of the Winger 
function has disappeared as W(0) = 0.075. 



dyne detectors are known for their high efficiency, the 
overall detection efficiency was about 75% in Wenger et 
at 's experiment 0] . An imperfect detector is equivalent 
to a perfect detector with a beam splitter in front pfij . 
where the transmittivity of the beam splitter is deter- 
mined by the detection efficiency rj. From ref. [3(J, we 
note that the characteristic function for the signal field 
passing through a beam splitter where the other input 
port is served by the vacuum is 



Substituting Eq. (|2*0j) into Eq. l|2*i)l. we find the char- 
acteristic function for the detected field. The Fourier 
transform of the characteristic function shows its largest 
negativity at the origin of the phase space and the value 
there is 

W im (0)oc-^=- ; 1 = (25) 



R(A-l) 



W 



R(B-l) 
1 



where v = T(A - l)r/ + 1 and w = T(B - l)r/ + 1. Under 
the assumption (A— 1)(B— 1) < 0, this becomes negative 
when the detection efficiency satisfies 

1 1 R 

V > ~2T(A- 1) ~ 2T{B- 1) ~ 4T' 

In particular, for a pure Gaussian incoming field, the 
condition becomes 

1 + T 



4T 



(26) 



(24) 



The RHS is smaller than unity (the detection efficiency 
X] < 1) only when T > 1/3. This is in good agreement 
with the perfect detection case. So, in order to see neg- 
ativity in the Wigner function, the beam splitter has 
to have a transmittivity larger than 1/3 first and then 
the detection efficiency has to satisfy the condition l|26|l . 
Wenger et al. employed a beam splitter with T w 0.88 
in which case the detection has to be larger than a mere 
53.4% to see negativity in the Wigner function. 

Another important factor which degrades the quantum 
effect of the photon subtracted Gaussian state in a real 
experiment is the modal purity factor |T^ . If the dark 
count rate of the photodetector employed to subtract a 
photon is non-negligible, the resulting state can be esti- 
mated in a mixture of photon subtracted squeezed state 
and squeezed state as 

tW(a) + (1 - 0W sg (a) (27) 

where W{a) is the Wigner function of the photon sub- 
tracted squeezed state, W sq (a) is the Wigner function of 
the squeezed state, and £ corresponds to the modal purity 
factor, which was 0.7, in Wenger et al.'s experiment 
The Wigner functions of the photon-subtracted Gaussian 
state have been plotted for a number of different cases in 
Fig. [21 It shows that the negativity of the Wigner func- 
tion disappears when both of the homodyne efficiency rj 
and the modal purity £ are considered taking relevant 
experimental values. We suggest that either the homo- 
dyne efficiency should be improved from 0.75 to 0.9 or the 
modal purity factor should be improved from 0.7 to 0.9 
to clearly observe the negativity of the Wigner function. 
In these cases, the minimum negativity will be —0.044 
and —0.073, respectively. 



IV. REMARKS 

In this paper, we are interested in the nonclassicality of 
a state produced by subtracting photons from a Gaussian 



6 



field. Subtracting a photon does not transform a classi- 
cal state into a nonclassical state whereas a nonclassical 
input remains nonclassical. This is in contrast to the 
case of adding a photon to a Gaussian field, in which 
case even a very chaotic field transforms into a nonclas- 
sical state [2a, I22I I23 . The non-Gaussian state obtained 
by subtracting a photon from a Gaussian field may show 
large negativity in its Wigner function. The condition 
to obtain the negativity is analyzed for a realistic case 
including the mixed-state input, threshold detection, in- 
efficient homodyne detection, and modal purity. The 
non-Gaussian state analyzed in this paper is compared 
with a coherent-state superposition which may be ex- 
tremely useful for fundamental and application reasons. 
The comparison shows the fidelity higher than 90% for 
the experimentally relevant situation. We compare our 
analysis with a recent experimental demonstration |l3| 



of the photon-subtracted Gaussian field and suggest that 
either the homodyne efficiency or the modal purity factor 
should be improved to around ~ 0.9 to clearly observe 
the negativity of the Wigner function. 

Note added After completion of the work, we were 
made aware of ref. 31] which considers nonlocality of 
a photon subtracted squeezed state. 
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